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Abstract. A thorough analysis of stability and perturbativity bounds is performed in several 
jQ^' versions of the two-Higgs doublet model, for both CP-conserving and spontaneously broken CP 

minima. LEP results further aid in establishing very strict constraints on the mass of the lighter 
Higgs particle. 
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^ ■ Despite the great successes of the Standard Model (SM) of particle physics, it leaves many 

i unanswered questions, such as the origin of matter-antimatter asymmetry; although the SM does 

contain a CP-violating parameter in the CKM matrix, and violates baryon number, it is generally 
accepted that it does not lead to baryogenesis sufficient to explain the observed asymmetry. One 
of the simplest extensions of the SM, which tries to solve this problem, is the two-Higgs doublet 
I model (2HDM) p!], wherein a second Higgs doublet is added to the theory. The spectrum of scalar 

particles becomes richer and, for some realisations of the model, spontaneous breaking of the CP 
symmetry is possible. The 2HDM presents some challenges, though: except in super symmetric 
models, the quartic interactions between the scalar doublets are not theoretically constrained, 
and increase substantially the number of free parameters. As a consequence the predictive 
power of the model is reduced. Any tool available to constrain the parameter space is thus 
of great interest. In this paper we will take a closer look at the requirements of stability 
and perturbativity of the model. Namely, we will analyse their impact on the several possible 
incarnations of the model; including versions involving the imposition of two types of global 
symmetries, which eliminate several unknown parameters. For various reasons, it may be of 
interest to break those symmetries softly, by the introduction of quadratic coefficients in the 
potential. The possible different vacua of the model - minima which spontaneously break CP 
or preserve it - require a separate stability and perturbativity analysis, which we will perform. 
This paper is organised as follows: in Section [2] we will briefly review the basic notions about 
the 2HDM scalar potential and the requirements of stability and perturbativity for a range 



*ferreira@cii.fc. ul.pt 
Mrtj@sune.amtp.liv.ac.uk 



of renormalisation scales. This will lead to the computation of the one-loop /3-functions of 
the model, (previously given in Ref. [3j). In sections [3] to [6] we will apply the stability and 
perturbativity bounds to the several realisations of the 2HDM: models with a discrete Z2 or 
global symmetries and their softly broken counterparts; within these, we will consider the 
possible cases of minima with spontaneously broken CP, or unbroken CP; and finally we will 
also consider the most general CP-conserving 2HDM potential. In all cases, we will endeavour 
to obtain bounds on the masses of the scalar particles, and use the latest experimental results on 
Higgs searches from LEP [2] to further constrain the potential's parameter space. Details of the 
/3- function calculation are given in Appendix [Aj following a simple and pedagogical approach 
which may be of interest for readers unfamiliar with it; and in Appendix[B]we make some remarks 
about the renormalisation group invariance of basis-invariant conditions on the couplings. 



2 The 2HDM potential 

The 2IIDM potential [1] involves two Higgs doublets with hypercharge Y = 1, $1 and ^2, and 
is invariant under the gauge symmetries of the standard model, SU{3)c x SU{2)w x {7(l)y. 
The most general potential one can build with these two doublets, following the conventions 
of [5], is given by 

V =mli^\^i + 77122^2*^2 - ("T'12^1<^'2 + h.C.) + 



^ + y ($5^2)' + A3 ($l$l) (^$2) + A4 ($1^>2) (^1^1] + 



|y (^1^2)' + [As + Ar {^l^, 



^>|$2 + ll.C. 



where the couplings {?ti^2 1 A5 , Ae , A7} are in general complex. In all, this potential has at 
most 14 real parameters. If one defines the CP transformation of the scalar fields as $1 — > <I>J, 
<I>2 *I*2' requires that the potential above preserves this symmetry, all the parameters 
become necessarily real. The number of free real parameters is thus reduced to 10, or in fact 9, 
with an appropriate choice of basis for the scalar doublets. Namely, a given linear combination 
of <I>i and ^2 will always diagonalise the quadratic terms in the fields in Eq. ([T]). For all that 
follows, we will consider that the potential does not break CP explicitly, and thus all parameters 
are taken as real. 

In general both $1 and $2 could have distinct Yukawa couplings to up-type quarks, down- 
type quarks and leptons. However, these generic Yukawa couplings would induce flavour- 
changing neutral currents (FCNC) in the theory. These have to be kept in check, either by 
imposing severe bounds on the size of the model's parameters or, more elegantly, by imposing 
symmetries upon it. Namely, a discrete Z2 symmetry [6J or a global U{1) [7] will prevent any 
FCNC from arising. This can be accomplished in several ways, but the choice we will make in 
this paper is to have only <I>i coupling to fermions. The results for the /3-functions of the model 
are easily generalised to other situations, by means of the techniques detailed in Appendix lAl 
In what follows we will, however, retain only the top quark Yukawa coupling, as the remaining 
ones will be too small to have any meaningful effect on the analysis we will perform. With 
this assumption we present in Appendix |A] the one loop /^-functions for the theory defined by 
Eq. ([1]). Expressions for the /3-functions for general models may be found in the literature [8], 
and the explicit expressions including the Xqj contributions were given in Ref. [3]. Even with 
FCNC-preventing symmetries imposed upon it, the 2HDM potential has a great number of free 
parameters - 7 or more of them - a fact which severely curtails its predictive power. Any tools 
which help in limiting this vast parameter space are thus welcome. One way to limit the values 
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of the quartic couplings in Eq. ([T]) is by observing that general values for the Aj do not guaran- 
tee that the potential is bounded from below (BFB). In fact, lest one requires that the quartic 
terms in Eq. ([T]) do not tend to minus infinity for any direction in field space, we will have no 
guarantees that the potential can have a stable minimum. For potentials where As = A7 = 0, 
Ivanov |9j has proven that the 2HDM potential is bounded from below if and only if the following 
conditions are obeyed: 

Ai > A3 > - VA1A2 

A2 > A3 + A4 - IA5I > - . (2) 

These conditions have been widely used in the literature and assumed to be only necessary ones, 
but they are also in fact sufficient. The work of [9] gives, in principle, all necessary and sufficient 
conditions to have the potential Eq. ([T|) bounded from below even in the case Ag / 0, A7 7^ 0, 
but the relations one could derive in this situation are extremely complicated, and not at all 
clear. See also [10]. In [16] necessary conditions involving Ag and A7 were derived, and we will 
use them in the following work: 

2|A6 + A7I < + A3 + A4 + As . (3) 

As explained above, the conditions Eqs. ([2l) and ([3]) ensure the stability of the tree scalar 
potential. To be sure of a viable vacuum, however, one must take into account the effect of 
radiative corrections, and the related fact that the Aj depend on the renormalisation scale 

Let us first review this important issue in the context of a theory with a single scalar field, 
the real scalar (/i^-model, with 



Vd = + (4) 



so that 

T/ = T/ , 4- - 

647r2 V /"^ 

Let us suppose for simplicity that > 0. For what values of (f) can we reliably calculate VI 
Suppose we have chosen a RG scale /i ~ m and that A is small on that scale so that perturbation 
theory in A is believable. Then evidently we can calculate V as cj) ^ Qhy simply retaining /i ~ m 
since the one loop correction is obviously small. Thus the origin remains a minimum, as was 
the case for the tree potential. But what about (j) » ml The one-loop correction now becomes 
large, because of the logarithm, so that one must improve on this perturbation expansion. RG 
improvement amounts, in fact, to exploiting the freedom to choose the renormalisation scale to 
take ji^ ~ V"J(0), or /i ~ for large 0. Then to a good approximation at large we will have 

V = ^0* (6) 

and this will be perturbatively believable as long as \{(t)) is small. Now in this simple model 
A becomes large at large scales, approaching a Landau pole, and so perturbation breaks down 
eventually in spite of our RG improvement. Thus we cannot say what form the potential takes 
at sufficiently large (p. 

In a more complicated theory there are two main issues to take into consideration. Firstly, 
if the potential depends on more than one scalar field, it is not immediately obvious in which 
directions in field space we will be able to describe the large-field potential, since we have 
only one scale at our disposal- Secondly, the behaviour of the Aj((/)) for large cf) may be quite 
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different from that in the simple i;^^-model. In particular, in the 2HDM the large size of the top 
quark Yukawa coupling, and the sign of its contribution to the /3-function of Ai (see Eq. (|44p in 
Appendix |A]) drives down the value of that quartic coupling as one goes up in renormalization 
scale. If the starting point of Ai is sufficiently small, Ai may become negative at a given high 
value of /i, which would mean that any minimum present for low renormalization scales would 
in fact be unstable - the potential would either be unbounded from below or develop a much 
deeper minimum at large (f)^. 

The approach we shall take is to simply assume that the stability conditions Eqs. ([2]), ([3]) must 
hold at all renormalisation scales fi up to the (putative) gauge unification scale Mu = 10^^ GeV. 
This will clearly be sufficient to produce a potential bounded from below. Requiring the stability 
of the scalar potential at all scales will thus, typically, impose lower bounds on the values of its 
quartic couplings. 

Another way of limiting the values of the Aj is by requiring that they remain small enough 
for perturbative believability at high scales. Hence, if the initial values of Aj are too large, their 
/3-functions will be positive and their renormalization scale evolution will drive them to ever 
higher values. Requiring that the Aj remain small at all scales will thus impose upper bounds 
on their values. How small should "small" be? Here we enter a somewhat arbitrary region, 
but requiring that all Aj remain less than 10 at all renormalization scales seems a reasonable 
requirement. 

We will therefore impose both stability and perturbativity bounds on the quartic parameters 
of the 2HDM at all scales between the weak scale Mz and Mfj. Such analyses have been made 
before, in many works: these ideas were applied to the SM [18], SUSY models [19] and also to a 
simple 2HDM [20] . In this work we are interested in studying the differences that the application 
of these bounds will have on the several possible two-Higgs doublet models, and on the several 
possible vacua therein possible. 



3 Model with Z2 symmetry 

One of the symmetries that rids the potential Eq. ([T|) of FCNC was first proposed by Glashow, 
Weinberg and Paschos [6j, and consists of a simple Z2 transformation in the fields: $1 — > — ^"1, 
$2 <^2- By carefully choosing similar transformations for the fermionic fields it is possible to 
eliminate the existence of FCNCs by having, for instance, only couple to the fermions. This 
symmetry simplifies Eq. ([1]), namely setting to zero several of the couplings: mi2 = Ag = Ay = 
0. Then the BFB conditions of Eq. ^ are, in this case, necessary and sufficient. 

Our procedure was as follows: we generated many thousands of combinations of quartic 
parameters of the 2HDM. The couplings were generated with magnitudes between 10^^ and 10, 
allowing different couplings to have different orders of magnitude and to be negative if allowed. 
Spontaneous symmetry breaking occurs when the doublets acquire vacuum expectation values 
such that 

(") • ■ m 

We thus generated values for the vevs {vi , ^2} such that + V2 = v"^, with v = 246/\/2 
GeV. With the vevs and all quartic couplings, it is simple to use the stationarity conditions 
of the model and determine the quadratic parameters mli and m22- At this point, we have a 
full set of parameters for the potential. By analysing the model's squared scalar mass matrices 

^In fact, in the type of theory we consider here, the latter is generally the case because the positive contribution 
of the gauge coupling contributions to their /3-functions causes Ai to recover to positive values at yet higher scales. 
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(expressions for which may be found, for instance, in [16j). we can ensure that each combination 
of parameter values under consideration is indeed a Normal minimum of the 2HDM. 

We then analysed the RG evolution of the quartic couplings for each "point" of parameter 
space and checked whether they obeyed the stability and triviality bounds described above, 
between Mz and Mu. In Fig. [1] we see the result of this procedure. In this plot we show the 




Figure 1: Results of scan of the 2HDM potential with Z2 symmetry, 
lightest CP-even Higgs mass versus the value of the Ai coupling. The colours are interpreted as: 

• The red (medium) points represent those combinations of 2HDM parameters for which the 
stability conditions of Eq. ^ were violated somewhere between Mz and Mu. 

• The blue (dark) points represent all of the parameter combinations which passed the 
stability conditions of Eq. (I2|), but for which no triviality conditions were set. 

• Finally, the green (lightest) points are a subset of the blue ones - those for which the 
triviality conditions are obeyed at all scales between Mz and Mjj H. 

As we see, the combination of stability and triviality conditions narrows the "allowed" range 
of Ai immensely - the only values which "survive" are in the interval 0.24 < Ai < 0.91. The 
remaining couplings are likewise constrained in similar intervals, of identical order of magnitude. 
This also limits high values for the Higgs scalar masses. In fact, if one analyses the full spectrum 
of scalar particles, one finds that after the stability and triviality requirements the masses are 

■^Notice that there are many blue points "between" the red ones, even if in the plot they are "covered" by the 
red points. This means that not all points with Higgs mass below ~ 200 GeV are excluded on the basis of the 
stability conditions, and that the green region is indeed a subset of the blue one. This is not a contradiction - this 
parameter space includes seven different parameters, and this plot is only varying two. We can therefore have a 
"rejected" point and an "accepted" one occupying the same place in the plot. 



bounded by (roughly): 



ruh < 102 GeV 
121 < ruH < 199 GeV 
niA < 163 GeV 
mH± < 160 GeV . (8) 

These results, specially those pertaining to the charged Higgs mass, are in agreement with 
previous works (see the last reference of [20]). Unless explicitly stated, no lower bounds were 
found for these masses, the exception being the heaviest neutral scalar H. These bounds do not 
preclude very low Higgs masses, then. In fact, current experimental data does not forbid light 
neutral 2HDM scalars. The best bounds on Higgs masses arise from the latest LEP results [4] 
and the analysis of associated production of a Z boson with the lightest 2HDM CP-even scalar, 
through the triple vertex Z Z h. In the 2HDM, the coupling associated with this vertex is equal 
to its SM value, multiplied by sin(a — /?), where tan/3 = V2/V1 and a is the usual mixing angle 
for the CP-even scalar mass matrix. If this coupling is small - meaning, if sin(a — f3) is small 
- then the mass of the scalar particle can be also small and still have escaped detection at 
LEP. This was explored recently for both SUSY models and one version of the 2HDM (see, for 
instance, [21]). In fact, the cross section for e~ — > Z h production in the 2HDM is related to 
that of the SM by 

cj2™^(e+e- ^ Zh) = sin2(a - /3) a^*^(e+e- ^ Zh) . (9) 

This relation is valid for any type of 2HDM model with a CP-conserving vacuum. The LEP 
results impose severe constraints on the size of the ratio (j'^^^M ^^SM considering the 

previous equation, translate as constraints on sin^(a— /3). In Fig.[2]we plot the value of sin^(a— /?) 
against the mass of the lightest Higgs boson, for the subset of parameter space which survived 
the stability and triviality bounds (the green (light) points from Fig. [1]). The red (continuous) 
line drawn in this plot corresponds to the experimental bound coming from the LEP searches 
Only the points below this line are allowed. We therefore see that the majority of the points 
which survived the triviality and stability analysis are already excluded on experimental grounds. 
Although it seems possible to generate high masses for the lightest Higgs particle, this plot clearly 
indicates that lower masses are preferred. Indeed, a rough upper bound of ~ 55 GeV can be 
established from these data. 

Nevertheless, caution must be urged. However large our sampling of the parameter space, 
it does not cover all regions of it. Also, these results are sensitive to the input top quark mass, 
which still has some uncertainty, according to the most recent Tevatron results |!23j: the CDF 
and DO combined value is Mt = 173.1 it 1.3 GeV. The physical mass corresponds to the pole 
of the propagator, and its relation to the Yukawa coupling ht and the vev vi is given, up to one 
loop, by 



Mr'' = hv, ^1 + (^4 - 3 In ^-^j ^JJ , (10) 

where we are only taking the most significant corrections, those from QCD. /x is the renormal- 
ization scale considered, and all quantities in the formula above are evaluated at that scale. The 
results presented thus far (and elsewhere in this paper) assume a top pole mass of 173 GeV. 
We verified what changes occur if we varied the top pole mass by 2 GeV in either direction (a 
conservative variation). The bounds shown in Eq. ^ that change by variation of Mf°^'^ are 
shown in Table [TJ The lower bound on the heaviest CP-even scalar is the one that changes the 
most. In fact, that lower bound correspond to small values of ^2, for which one of the masses, h 
or is essentially proportional to Ai. The lower bounds for mn presented in the table above 



LEP constraints for the potential with a symmetry 
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Figure 2: sin^(a — f3) versus the mass of the hghtest Higgs boson, for a potential with a Z2 
symmetry. 
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Table 1: Bounds on scalar masses in function of the value of the top quark pole mass. 

correspond to the lower allowed values for Ai, which obviously change when the top pole mass 
is varied. The uncertainty on the top pole mass is thus relevant, and needs to be factored in 
evaluating whatever bounds we will present in this work. 

Still, the results shown in Fig. [2] clearly indicate that the 2HDM with a Z2 symmetry is 
already severely constrained by the simultaneous requirements of stability, triviality and com- 
pliance with existing experimental results. The pole mass dependence is of the order of ~ 5 
GeV around the central values at most, and will not drastically change those conclusions. 

3.1 The case V2 = 

For the 2HDM with an unbroken Z2 symmetry, the minimisation conditions admit a different 
type of solution than the one we have been considering: to wit, a vacuum where one of the fields 
$ has a vanishing expectation value. These models were first proposed in [IT] and have been 
studied before, in many different contexts. For instance, in ref. [12], one such model was used 
to show that it was possible to have neutrino mixing even without massive neutrinos. In [13] 
the model was used to explain low neutrino masses as a loop effect. In general these models are 



the basis of the so-called "inert Higgs" theories [14J, which have excellent scalar candidates for 
dark matter [15j . 

Since in our models only <I>i couples to the fermions, we should therefore study the case 
where V2 = 0, and verify what changes occur in the bounds we have deduced 0. The expressions 
for the squared scalar masses in these models are extremely simple, namely 

m]i± = ml^ + Xsv^ Ml = m^^ , (11) 

where m| and m| are the two CP-even scalar masses, the lightest of which will be h, the 
heaviest H (depending on the parameters, though, we cannot a priori guarantee which of hi 
and /i2 is the lightest state). In the model of the previous section, the minimisation conditions 
ensured a strong bond between the values of the squared parameters, mfi and m22, the scale v'^ 
and the values of the quartic parameters Aj. With V2 = 0, though, despite the fact that mli 
is fixed such that mf^ = — Aif^, the parameter 777-22 is unconstrained and as such can be as 
large as one wishes. Thus, we expect the upper bounds on most of the masses written above to 
be much larger than before. In fact, once the stability and triviality analysis is concluded, we 
obtain 

mh < 235 GeV 

mn > 120 GeV . (12) 

The only upper bound that remains is that on ruh, certainly due to the hi state, which is directly 
tied to the severely constrained Ai coupling. Because the CP-even 2x2 mass matrix is diagonal, 
the mixing angle a has only two possible values: and tt/2. As such (and since in this model 
the angle /? is equal to zero), the LEP results have no impact on the model: 

• If a = 0, the coupling of h to the Z boson vanishes. This occurs for a large range of 
masses, from very low Higgs masses to high ones. In any case, the LEP data do not 
provide any constraints, since this lightest scalar does not couple to the Z and as such 
could not have been observed at LEP. 

• If a = 7t/2, the coupling of /i to Z is identical to that of the SM. However, this case only 
occurs for values of the Higgs mass larger than 116 GeV. The LEP constraints, as can be 
observed from fig. [21 are only valid for Higss masses inferior to about 110 GeV. As such, 
the case a = 7r/2 is also not constrained by the LEP data. 



4 Model with softly broken Z2 symmetry 

If one adds to the Z2 potential a term of the form 771^2*^1*^2 + h.c. the discrete symmetry is 
softly broken. However, no FCNCs arise from this soft breaking. The main reason why doing 
this should be of interest is quite simple - with this soft breaking term the potential can now 
have two types of interesting minima: (a) "normal" ones, which preserve CP, and for which the 
doublet's vacuum expectation values, as before, have the form of Eq. ([7j); and (b), minima for 
which CP is spontaneously broken, the doublets developing vevs of the form 

<*fi>= , <*f2>= f °. ) ■ (13) 



'We thank the referee for bringing this point to our attention. 



In order to have minima of phenomenological interest, the CP vevs must obey vf + V2 + v'^ = 
GeV^. It has recently been proved [9] that these different minima cannot co-exist. The 
parameters of the potential either allow the normal minimum or a CP breaking one i Thus, 
the parameter sets for each of these cases have to be different, and generated separately. The 
different vevs, however, do not directly affect the RGE running of the quartic couplings. Notice 
also that the existence of another free parameter {vn\2l the potential will certainly change 
the scalar masses. 

4.1 Normal minima 

For this case, the parameters of the potential were chosen such that the global minimum of the 
potential preserves CP. The range of allowed values for the Aj was the same as before, and the 
soft breaking parameter is chosen so that all three mass squared parameters are of the order 
^ and so that \'rn\2\ < Wv'^. Once again, thousands of different parameter sets were generated, 
and then RGE analysed between the weak and unification scales. Again we found that only a 
narrow window of values of Ai survived the imposition of stability and triviality bounds. The 
results differ from those of the unbroken Z2 model, the bounds found for the masses being given 
by 

ruh < 187 GeV 

121 < niH < (no upper limit) GeV . (14) 

As explained earlier, the quadratic parameters and 17122 obtained from the stationarity 
conditions of the potential, once all other parameters have been generated. For the unbroken 
Z2 potential, those included all the quartic couplings (limited by the triviality and stability 
requirements) and the vevs (limited by the requirement of their squared sum be equal to (246/2)^ 
GeV'^). Now, however, we have the extra quadratic parameter in the potential, which is not 
restricted and can make the masses larger than they can be in the unbroken Z2 model. This 
justifies the lessening of the bounds we discover. 

In Fig. ([3]) we again plot the effect of the LEP bounds; once again the allowed region is 
below (to the right of) the line. As we see, the soft-breaking term allows us to easily evade the 
experimental constraints. 

4.2 CP breaking minima 

As was shown in Ref. [16], the quartic parameters of the 2IIDM potential need to obey a 
specific condition so that there is a CP-breaking minimunj^l. Thus, if the quartic parameters 
are necessarily different from those which generate Normal minima, one expects different results 
stemming from the stability and triviality bounds. 

Notice, now, that there is no distinction between CP-even and CP-odd scalars, since CP 
is spontaneously broken - the complex vevs shown in Eq. (jl3p will cause a mixing between 
all neutral components of the doublets. We rename the neutral scalars as hi, h2 and /13, in 
decreasing order of masses. After repeating the stability and triviality RG analysis, we found 
the following bounds for these scalar masses: 

131 < ruh^ < 203 GeV 

^Charge breaking vacua might also occur but, as shown in [16j . whenever a normal minimum exists in the 
2HDM, the global minimum of the potential is normal and thus safe against charge breaking . Likewise, if a CP 
minimum exists, any charge breaking stationary point will necessarily be a saddle point. 

^Namely, the quartic parameters must be such that the matrix Bcp defined in [TS] be positive definite, and 
so that A4 < A5. 
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LEP constraints for the softly broken potential. Normal minimum 
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Figure 3: sin^(a — /?) versus the mass of the hghtest Higgs boson, for the normal minimum of a 
potential with a softly broken symmetry. 

6 < m,,2 < 160 GeV 
< 84 GeV 

m^i < 158 GeV . (15) 

Notice the very low upper bound on the lightest neutral scalar, and indeed in all of the scalars. 
Clearly the requirement of a CP-breaking vacuum "chooses" , from within the parameter space 
which survives the stability and triviality bounds, scalars with lower masses, which is quite 
different from what we saw occurring for this same model, for Normal minima. One might 
speculate that this phenomenon is somehow related to the Georgi-Pais theorem [26] according 
to which spontaneous breaking of CP via radiative corrections is always accompanied by scalars 
which are massless in the tree approximation (of course here we are considering tree-level CP 
breaking). 

The LEP results also have a substantial impact on the parameter space for these minima. 
But to apply them, we must first compute the coupling between the Z boson and the lightest 
Higgs scalar for a CP-breaking minimum. As was explained earlier, the complex vevs of Eq. (|13p 
cause a mixing between CP-even and CP-odd scalars, so that the masses of the neutral scalars 
/ii, /i2 and /i3 are the eigenvalues of a 4 x 4 matrix (the fourth eigenvalue is zero, corresponding 
to the Z would-be Goldstone mode G*^). Consequently there is now no single angle a which 
characterises the diagonalisation of this matrix, and thus the quantity sin(a — /3), which was the 
ratio between the Z Zh coupling in the 2HDM and the SM, is no longer defined for CP-breaking 
minima. For the determination of the Z Zh coupling, the relevant term in the Lagrangian stems 

^We thank the referee for a comment on this point. 
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from the kinetic terms for the doublets, namely 



i=l 



\ sec^ 



Let us define the neutral component fields of the doublets as 



(16) 



(17) 



The relationship between these fields and the mass-eigenstates /ii, /i2, /i3 and is given by a 
4x4 unitary matrix Ajj, such that 



(18) 



Then, the terms in Eq. (jl6p which are directly proportional to the lightest Higgs, /13, will be 
given by 
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so that the ratio of the Z Z h coupling in a CP minimum and that of the SM is given by 

gzzh = - (Aisvi + A23V2 + A43V3) 

V 

with V = 246/\/2 GeV. Therefore, the generalisation of Eq. ([9]) is thus 



2HDM . 
^CP minimum'' 
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It is now a simple task to calculate the mass matrix for the neutral scalars, diagonalise it and 
obtain the matrix A and thus compute the value of gzzh for the points of parameter space which 
survived the stability and triviality bounds, comparing them with the LEP results. As we see 
from Fig. ([5]), the LEP results exclude a significant portion of the parameter space. (As before 
the allowed region is below the line). In fact, according to this plot, the highest value allowed 
for the lightest Higgs mass would be around 65 GeV. Comparing this plot to figs. [2] and [3l we 
see that the bounds we are imposing have very different consequences, depending on the model, 
or type of minimum, considered. 



5 Model with a softly broken ViV) symmetry 

Another symmetry to eliminate FCNC in the original potential Eq. ([T]) is a simple global C/(l) 
transformation in the fields (accompanied by suitable fermion transformations) of the form 
$1 — > $1, $2 — *■ e*"^'2. This symmetry simplifies Eq. ([T]), again setting to zero several of 
its couplings, m\2 = A5 = Ar = A7 = 0. This symmetry is however too strong, in that 
it produces a zero mass axion To prevent that from happening one usually softly breaks 
this global C/(l) by re-introducing the 771^2 term in the potential. Except for the fact that the 
remaining A couplings are unrelated to one another, the resulting Higgs potential is similar to 
the MSSM one. The only types of minima this potential possesses are normal ones; CP breaking 
is impossible here (see also footnoted]). 



Unless one considers a vacuum with V2 = 0, as in section [3Tll with analogous consequences. 
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Figure 4: g^zh versus the mass of the hghtest Higgs boson, for the CP breaking minimum of a 
potential with a softly broken Z2 symmetry. 

Once again we generated thousands of parameter sets corresponding to normal minima 
and used the model's /3-functions to verify whether the stability and triviality conditions were 
satisfied all the way up to Mjj- In terms of that analysis, the only difference with the softly- 
broken Z2 model studied in section 14.11 is the fact that for this model we must have A5 = 0, 
even after the U{1) symmetry has been softly broken. As before, only a very narrow range of 
vales of Ai survived the stability and triviality requirements. The bounds found for the scalar 
masses are now: 

rrih < 187 GeV 

121 < niH < (no upper hmit) GeV . (22) 

The LEP constraints do not affect the available parameter space as much as in the model with 
unbroken Z2 symmetry, or for CP minima in the softly broken Z2 model; the results are shown 
in Fig. ([5]) . The results are quite similar to those obtained in section 14.11 for normal minima in 
the softly-broken Z2 model, no doubt due to the presence of in both cases of the soft-breaking 
parameter m|2- 

6 The full CP conserving potential 

What about the full potential of Eq. ([T])? What do the triviality and stability bounds tell us 
about it? To perform this analysis we need the /3-functions for this model, in terms of the "new" 
couplings Ag and A7. These are to be found in Appendix El As we have mentioned, in this case 
FCNC will, in general, occur. 

Due to the presence of the new quartic couplings, the stability analysis needs to take into 
account Eq. ([3]), which involves Ag and A7. And once more, as with the softly broken Z2 model. 
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Figure 5: sin^(a — /?) versus the mass of the Hghtest Higgs boson, for the CP breaking minimum 
of a potential with a softly broken U{1) symmetry. 

this potential can have minima with spontaneous CP breaking and Normal minima, though the 
same set of parameters cannot produce two such minima in coexistence. We thus obtain the 
following bounds for each type of minima: for Normal minima, the only bounds found were 

mh < 187 GeV 

121 < niH < (no upper hmit) GeV , (23) 

and for CP minima, 

131 < mh^ < 205 GeV 
6 < mh2< 166 GeV 
m/,3 < 84 GeV 

mH± < 158 GeV . (24) 

As we can see, there are no great differences with the results obtained for the Normal minima 
of the softly broken Z2 potential, or for the U{1) model. The new parameters do not affect 
the bounds on the masses, nor do they change the qualitative difference between the bounds 
obtained for each type of minima: that the requirements of stability and triviality tend to "pick" 
lower scalar masses for CP minima than they do for Normal minima. Also, notice that the case 
V2 = 0, discussed in section [STTl is only possible for the Z2 symmetric potential (in fact solutions 
with V2 = become possible as long as = Ae = 0, but without the full Z2 symmetry these 
conditions would not be preserved by renormalisation) . It corresponds to a completely different 
type of vacuum (one which preserves the Z2 symmetry) for which the bounds we found for that 
case (eq. (fT2]) ) are unaffected. 

The LEP constraints are again similar to those already shown, and we present them, for 
both minima, in Fig. [6l Once again, much of the available parameter space is excluded by the 
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Figure 6: Squared ZZh coupling versus the mass of the lightest Higgs boson, for the full CP 
conserving potential. The green (light) points concern the Normal minima, the blue (dark) ones 
the minima which spontaneously break CP. 

LEP data for the case of the CP minima. The lightest Higgs scalar, in that case, would have a 
rough upper bound of ~ 80 GeV. 

7 Conclusions 

We have performed a thorough analysis of the impact that the demands of stability and triviality 
have on the scalar masses of the 2HDM. We considered several possible incarnations of this model 
- models with Z2 or U{1) symmetries, with those symmetries softly broken or simply without 
them, and the different neutral vacua allowed in those theories. At the same time, we studied 
the impact of the LEP results on production of a light Higgs scalar on the parameter space of 
the model. Our results may be summarised as follows: 

• The 2HDM potential with a Z2 symmetry is very strongly constrained. Combining both 
theoretical and experimental constraints, the mass of the lightest neutral scalar should be 
less than about 55 GeV. 

• The LEP restrictions are easily avoided in Normal minima of models with softly broken Z2 
or U{1) symmetries, or in the full CP-conserving potential. In those models, the lightest 
CP-even neutral scalar is bound to be smaller than about 190 GeV. The heaviest CP-even 
neutral scalar is bound to be larger than about 120 GeV. 

• The minima with spontaneous CP breaking which may occur in these models are heavily 
constrained. For these minima, the stability and triviality bounds affect the parameter 
space of these models in different ways than what occurs for the Normal minima. Those 



bounds constrain very tightly all scalar masses, and the LEP results have an extremely 
strong impact on the surviving parameter space, eliminating most of it. 

In order to perform this analysis, we needed the /3- functions for the parameters Ag and Ay, 
which are given in Appendix |A] and were given previously in [3J. These /?- functions allowed 
us to verify the validity of necessary conditions involving Ag and A7 between the weak and 
renormalization scales. Their usefulness, however, is not restricted to the studies presented 
here. In [27], for instance, basis-invariant conditions which ensure greater symmetry of the 
2HDM potential were obtained. Several of those conditions involve the parameters Ag and Ay. 
To verify if these conditions are valid at all renormalisation scales, it will be necessary to employ 
the /3-functions. We perform this analysis in Appendix [Bl 

The parameters Ag and Ay are often omitted because in their presence there is no global 
symmetry which we can use to automatically prevent flavour changing neutral currents. We 
have avoided this issue by assuming that quarks couple to $1 only, and anyway retaining only 
the top quark Yukawa coupling. Of course the absence of a symmetry to enforce this means that 
to pursue the question of FCNCs we would need to consider the effect of radiative corrections. In 
any event, we found the effect of Agj on our analysis to be limited - the analysis of the stability 
and triviality bounds for the full CP-conserving potential did not produce significantly different 
results from those obtained for the softly broken Z2 or U{1) models. That in itself, however, 
is an interesting result: if for whatever reason one wishes to work with the full CP-conserving 
2HDM model, one can do so with the certainty that the Ag and Ay parameters will not spoil the 
restrictive bounds obtained in simpler models. 

Finally, the bounds we obtained here considered that the stability and triviality conditions 
held at all scales between Mz and Mjj. More to the point, this procedure has the underlying 
assumption that the two-Higgs doublet model constitutes the whole of physics up to the gauge 
unification scale. That philosophy, however, can be readily inverted. A possible scenario is that 
within the next few years several scalar particles are discovered at the LHC, and that their 
properties conform to the 2IIDM. However, suppose their masses completely violate all bounds 
presented here. This would, of course, suggest the existence of more new physics beyond the 
2HDM but below the gauge unification scale, to justify the breaking of the stability and triviality 
bounds. The simplest example of such new physics would be the existence of a heavy fourth 
family of fermions, the presence of which would significantly change the form of the /3-functions 
of the model, by virtue of the (necessarily large) associated Yukawa couplings. 

Acknowledgements: This work was supported in part by the Portuguese Fundagdo para 
a Ciencia e a Tecnologia (FCT) under contract PTDC/FIS/70156/2006. DRTJ also thanks the 
CERN Theory Division for hospitality and financial support. 

A The one loop /^-functions 

In this appendix we describe the calculation of the one loop /3-functions for the general two 
Higgs scalar model defined by the potential given in Eq. ([1]). The calculation is straightforward 
by normal diagrammatic methods, and has already been presented in Ref. |3j. Here we describe 
an alternative algebraic calculation based on the renormalisation group equation satisfied by the 
effective potential, which we hope may be of some pedagogic interest, as well as providing a check 
on the previous calculation; in which we indeed thereby identify one fairly obvious typo. This 
procedure uses the RG invariance of the one loop effective potential, where all the /3-functions 
may be computed from STrM^, where is the mass matrix of the fields, including arbitrary 
vevs for the scalars. The basic procedure is explained in section 6 of reference |24j . 
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Up to one loop, the effective potential Vesifj)) for any theory is given by: 

V,s{^)=V{^) + Vi{^) + --- (25) 
where V{(j)) is the tree potential given in our case by Eq. ([1]), and 



yi(^) = ^STrM^ln^ , (26) 
4 



where the mass matrix includes contributions to all scalar, fermion and vector boson mass 
matrices with arbitrary background values of all scalar fields cj), k = (IGvr^)"^ and STr is the 
usual spin-weighted trace. 



In the Landau gauge, V^fj obeys the following RG equation: 

d ^ „ d , . d 



^cfr = (27) 



where the Aj include all mass parameters and coupling constants, and 7 is the matrix of anoma- 
lous dimensions of the scalar fields. 

It follows that 

V^^)V = -^i^Vi = -STrM^ , (28) 
Ofi 2 

where 

^'"'-Erj|-(^V">|+c.o.). (2.) 

By comparing coefficients of the various ip'^ terms on the two sides of Eq. (j28p we can, if we know 
7, determine all the one-loop /3 functions. From now on we write P^^^^ = K,(3i and 7^"^^ = K7 to 
avoid writing factors of k. 

Let us first consider the simplified case when we set gauge and Yukawa couplings to zero. 
Then we can write as follows: 

M'=i ^ n 1 (30) 

where 



A = ^-^gf (31) 
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and i,j = 1,2 are SU{2) indices, and to control the profusion of indices we have put = cp 
and ^2 = i- 
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It is then straightforward to write down the matrices A, B, C. Thus, for example, 



(35) 



Since we are neglecting gauge and Yukawa terms we have no contribution to 7, we have 
from Eq. ([25]) that 



- STrM^ 
2 



^ (Tryl^ + 2Tr5C + TtD'^) 



+ 



13m 
2 

2 



+ 



+ h.c. 



(36) 



It is straightforward algebra to obtain from Eq. (I36p that 



/3ai 
Px, 
Pxs 
Px, 

Pxs 

Px, 

Pxr 



12X1 + + 4A3A4 + 2A| + 2A^ + 24A^ 
12A^ + 4A| + 4A3A4 + 2A^ + 2A^ + 24Af 

(Ai + A2)(6A3 + 2A4) + 4A^ + 2A^ + 2A^ + 4A^ + IGAeAy + 4A^ 
2(Ai + A2)A4 + 8A3A4 + 4Al + 8X1 + lOA^ + 4A6A7 + lOA? 
2(Ai + A2)A5 + 8A3A5 + I2A4A5 + lOA^ + 4A6A7 + lOA? 
I2A1A6 + 6A3(A6 + A7) + 8A4A6 + 4A4A7 + IOA5A6 + 2A5A7 
I2A2A7 + 6A3(A6 + A7) + 4A4A6 + 8A4A7 + 2A5A6 + IOA5A7. 



(37) 



We have verified the above results by a standard Feynman diagram calculation. Moreover they 
are in full agreement with the results presented in Appendix A of Ref. pj, except for a typo in 
the result for there; the contribution 12 Ag there should read 12 A7. (Note that there is an 
overall difference of a factor of 2 between the definitions of all the /3-functions). 

The contributions to STrM^ of 0{h'^) and 0{g'^,g'^g''^,g'^) are easily calculated. The 0{g'^,g'^g''^ 
terms come from the gauge boson mass matrix: 



{r^T^} 



+ 
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+ 







(38) 



Then using the identities 
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(39) 



one easily shows that 

5TrM^ 



+ 



^2/2 



)2 + (ete)2 + 20W^^ 



(40) 



where the STr has contributed a spin factor of 3. 

The 0{hf) contributions come from the top mass matrix: 



Ml 



(41) 



so that 



STrM/ 



-12ht 



(42) 



where the —12 consists of a colour factor of 3 and a spin factor of —4. 

The remaining contributions of 0(Aj/i^, Xig^, Xig'^) come from the anomalous dimension term 
in Eq. ()28p . the anomalous dimensions in the Landau gauge being 
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(43) 



Armed with these results one easily shows from Eqs. (j40p .()42 p . (j43|) that Eq. (j37p receives the 
following additional contributions: 



Px, 

Px, 
Px, 
Px, 

Pxr 



Px, + ^(35' + + 2g^g'^) - 3X^{3g^ + g'' - Ah',) - UK 
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Px, - SXriSg' + g' 



.12 



hi 



We also require the one-loop /3-functions for g,g',g3 and ht, which are given by 



(44) 
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B Renormalization group invariance and the exceptional region 
of parameter space 

As mentioned in section [2l one can impose various symmetries on the 2HDM in order to obtain 
interesting physical consequences. It was recently proved [28J that the 2HDM potential can in 
fact only possess six distinct symmetries (including both discrete and continuous symmetries). 
This statement, however, hides a problem: since physical predictions cannot depend on the basis 
chosen for the Higgs doublets, the form of the potential - its specific combination of parameters 
and respective values - is not uniquely determined. For instance, we discussed the 2HDM 
potential with a Z2 symmetry; nevertheless, a potential with a permutation symmetry between 
$1 and $2 - that is, a 2HDM potential invariant under the transformation $1 <-> $2 - has 
exactly the same physical predictions. The reason is that both potentials, and both symmetries, 
are related by a basis change on the scalar doublets (see, for instance, [29]). 

The question then arises, how does one know whether a potential has a given symmetry, 
since that symmetry can appear in an endless number of ways in different bases? The answer 
is, one builds basis-invariant quantities, the values of which reveal which, if any, symmetries the 
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potential has. There has been considerable attention to developing techniques to build basis 
invariants [30], first as a means of detecting CP violation and more recently to detect other 
types of symmetries jHH [27], such as U{1) or Z2. In fact, the authors of Ref. [31] built a set of 
basis-invariant quantities which can distinguish between these two symmetries. However, that 
method failed to identify the presence of a continuous symmetry for a given combination of 
parameters, the so-called exceptional region of parameter space (ERPS): namely, 

22 2 n 

"^22= "^11) = 0, 

A2 = Ai, A7 = -A6 . (46) 

Of particular relevance is the fact that the ERPS is not a zero-measure set of parameters, but 
is itself instead attained through the imposition, on the potential, of a given symmetry [31]. In 
fact, several possible symmetries lead into the ERPS: either combinations of discrete symmetries 
or generalised CP symmetries (generalised in the sense that they do not satisfy CP^ = 1). In 
fact, in [27J two generalised CP symmetries were identified, which lead into the ERPS: one was a 
discrete CP symmetry (dubbed CP2 in that reference), the other a continuous one (dubbed CP3). 
The problem of identifying the presence of a continuous symmetry was solved in Ref. [27] , with 
a new basis invariant quantity D, which is written, in the ERPS, in a basis where all potential 
parameters are real, as 

D = -^ [A5(Ai - A3 - A4 + A5) - 2XlY [(^1 - A3 - A4 - X5f + 16A2] . (47) 

As shown in [27j , once in the ERPS the basis-invariant condition D = indicates the presence 
of a continuous U{1) symmetry. If Ag = 0, then D = gives the following conditions: 

A5 = 0, A5 = ±(Ai-A3-A4). (48) 

If Xq 7^ 0, then D = corresponds to 

2A^ = A5(Ai-A3-A4 + A5). (49) 

In Ref. [27j the RG invariance of the D = condition itself was demonstrated, through a 
calculation which managed to avoid using the explicit form of the Aj /?- functions. Using the 
explicit form of the /^-functions of Appendix [Aj it is simple to demonstrate that: 

• The conditions on the quartic couplings that define the ERPS are RG invariant. Explicitly, 
we see that: 

— If Ai = A2 and Xq = — A7, then we will have (3\-^ = I3\^ and /^Ab = —Pxt 

• Each of the conditions in Eqs. (I48p and ()49p are RG invariant, if we are in the ERPS (that 
is, with A2 = Ai and Ay = —Xq). Namely, 

— If A5 = Ae = 0, then /Ja^ = Px, = 0; 

— If Ae = and A5 = ±(Ai — A3 — A4), then we have j3\^ = it {j3\^ — (3x^ — (3x^) and 
Px, = 0; 

— Finally, and much in the same manner, if = 2Ag — A5(Ai — A3 — A4 + A5) = then 
we also have f3^ = 0. 

However, there is a detail which must be mentioned: the above is true if one sets the Yukawa 
coupling ht equal to zero - that is, if the theory does not couple to fermions (the gauge coupling 
contributions from Eq. (I44p may be included, however). In fact, it is extremely difficult, if not 
impossible, to couple the two doublets to fermions in a phenomenological acceptable manner, 
if the symmetries which lead to the ERPS are in place. In Ref. [32j, for instance, the authors 
managed to couple the fermion sector and the scalar one in the presence of a CP2 symmetry, 
but those couplings implied masslessness for the two first generations. 
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